Universal phase diagram of disordered bosons from a doped quantum magnet 
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A quantum particle cannot in general diffuse through a disordered medium because of its wavelike 
nature, but interacting particles can escape localization by collectively percolating through the 
system. For bosonic particles this phenomenon corresponds to a quantum transition from a localized 
insulator phase - the Bose glass - to a superfluid phase, in which particles condense into an extended 
state. Here, we construct a universal phase diagram of disordered bosons in doped quantum magnets 
for which bosonic quasi-particles are represented by magnetized states (spin triplets) of the quantum 
spins, condensing into a magnetically ordered state. The appearance of a Bose glass leads to strong 
measurable signatures in the onset of superfluidity of the spin-triplet gas, exhibiting a complex 
crossover from low-temperature quantum percolation to a conventional condensation transition at 
intermediate temperatures. 

PACS numbers: 03.75.Lm, 71.23.Ft, 68.65.Cd, 72.15. Rn 



In continuous quantum phase transitions [1] quantum 
fluctuations drive a complete change of symmetry in 
the state of the system via the development of long- 
wavelength modes at arbitrarily low energies. This oc- 
curs in the case of strongly interacting bosons on a lat- 
tice, modeling widely different physical systems such as 
cold atoms loaded in optical lattices [2J, Cooper pairs in 
Josephson-junction arrays [5] , or magnetic quasi-particles 
in gapped quantum magnets H). I n tne presence of strong 
repulsion among the particles, a commensurate filling 
of the lattice stabilizes a gapped Mott insulating phase, 
which has all the symmetries of the original Hamiltonian. 
In contrast, doping the system away from commensurate 
filling drives it into a Bose-condensed phase [5], charac- 
terized by superfluidiy and a gap less phonon-like mode. 

The presence of disorder in the system can deeply affect 
its quantum critical behavior. In fact, disorder can local- 
ize the long-wavelength modes developed at the quantum 
phase transition of the clean system, and thus introduce 
a novel phase, with gaplcss localized excitations but in 
the absence of any long-range order. Indeed, it has been 
shown that lattice disorder for interacting bosons can lo- 
calize the gas of doped particles/holes, hindering their 
condensation and thus giving rise to a new insulating 
gapless phase with incommensurate filling, known as a 
Bose glass [5j |6]. Moreover, when the density of doped 
particles/holes exceeds a critical threshold, their mutual 
interactions destabilize energetically the localized state 
in favor of a delocalized superfluid state. It is clear that 
the nature of this transition is completely different from 
that of the Mott-insulator/superfluid transition in the 
absence of randomness. 

Despite its extreme richness and generality, the phase 
diagram of so-called dirty-boson systems has so far been 
extremely elusive, both to experiments, due to the diffi- 
culty of accessing the quantum critical regime in a con- 



trolled disordered environment [3 E] ; and to theory, 
due to the absence of well-controlled approximations in 
presence of strong interactions and strong disorder [TO] . 
Here we propose the use of doped quantum magnets to 
fully investigate the physics of dirty bosons, and we per- 
form numerically exact simulations on a realistic model 
Hamiltonian to construct the universal phase diagram 
of the system. Our model of interest is the magnetic 
Hamiltonian of NiCl 2 -4SC(NH 2 ) 2 (DTN) [THUS], which 
consists of coupled S — 1 chains with strong single-ion 
anisotropy: 

H = J c ^ (itjSi ■ Sj + J a b 2^ ei£m ^i ' Sm 

(ij)c (lm)ab 

+ dY,£*(SD 2 ~9VbHJ2^. (1) 
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where J c > is the dominant antiferromagnetic coupling 
for bonds (ij) c along the crystallographic c-axis, J a b > 
is the coupling for bonds (lm) a i, in the transverse a, b di- 
rections of the cubic lattice [T3] , and D > is the single- 
ion anisotropy. We adopt the experimental parameters 
HH EE2] Jc = 2.2 K , J ab = 0.18 K, D = 8.9 K. In the 
following we will use reduced temperatures t = UbT/Jc 
and reduced magnetic fields h — g^sH/Jc- Moreover 
we introduce randomness in the system in the form of 
site dilution via the variables q taking values and 1 
with probability x and 1 — x respectively, where x is the 
dilution concentration. Site dilution corresponds experi- 
mentally to Mg 2+ or Cd 2+ doping of the Ni 2+ ions. 

A standard spin-to-boson mapping [15) allows to recast 
the above Hamiltonian in the form of bosonic operators 
hi, b\, with the local Hilbert space limited to the Fock 
states |n) = |0), |1) and |2) corresponding to the states 
\ms)= | — 1), |0) and |1) of a S = 1 spin. Upon this 
transformation the J couplings take the role of nearest- 
neighbor hopping amplitudes and repulsion potentials, 
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FIG. 1: Phase diagram of coupled anisotropic S = 1 chains 
in a magnetic field with 15% dilution. We plot the critical 
temperature t c for magnetic BEC and the spin gap S = A/J c 
|18| as a function of the applied magnetic field for zero doping 
(x = 0) and for 15% doping (x — 0.15). The spin gap is inde- 
pendent of the doping concentration. The inset shows a zoom 
on the region close to the lower critical field h c i , marking the 
appearance of the Bose glass. 

while the single-ion anisotropy plays the role of on-site 
repulsion for the bosons, and finally the magnetic field 
acts as the chemical potential. 

Due to the large anisotropy D, the ground state of 
the system in zero field is well approximated by a col- 
lection of singlets \ms = 0), corresponding to a gapped 
Mott insulator of bosons with n = 1 particle per site. 
In the case of the pure system applying a mag- 

netic field above a critical value — 1.65(1) allows to 
close the gap in the spectrum and to increase the particle 
density above unit filling, corresponding to the appear- 
ance of a finite magnetization along the field direction. 
Upon increasing the filling, the system develops a super- 
fluid and condensate fraction, which translates into the 
spin stiffness of the magnetic system and into a peak 
at q = (-7T, 7T, 7r) in the transverse static structure factor, 

S±(g)=^ i exp[ig.(r i -r J 0]<^ ( " ) ^ (B V[(l-»)i 3 ],»- 
spectively. The condensation order parameter is directly 

related to the spontaneous magnetization appearing in 
the xy plane [16], (h) « (Sf - iS?)/(y/S - Sf). Con- 
densation (in the form of magnetic order) survives at 
finite temperature up to a critical temperature t c , which 
scales with the applied field as t c ~ (h — h 1 ^)^ , where 
(pa = 2/3 as predicted by mean-field theory for the con- 
densation of a diluted Bose gas jTTj. The experimental 
verification of this prediction [16] is considered as the 
strongest evidence for magnetic Bose-Einstein condensa- 
tion (BEC) . Increasing the field even further reduces the 
state of the system to that of a dilute gas of singly occu- 
pied sites (excess triplet holes), associated with \m s = 0) 
spin states. A high field eventually leads to full polar- 
ization of the spins at the critical value = 8.69, cor- 
responding to a second Mott insulating state of triplets 
(namely the vacuum of triplet holes), and to a vanishing 
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FIG. 2: Local magnetization profile (corresponding to the ex- 
cess spin-triplet bosons) of coupled S = 1 chains, Eq. |l| on 
a 16 3 lattice with 15% dilution, obtained via quantum Monte 
Carlo simulations. Four representative points in the h - t 
plane are shown [(h, t) coordinates are indicated below each 
figure], for the following bosonic phases: (1) a Bose glass with 
localized incoherent excess spin quasiparticles; (2) an inhomo- 
geneous BEC phase, with long-range coherence due to weak 
links between localizes states; (3) a normal Bose gas, with 
incoherent spin quasiparticles which are spread over a perco- 
lating region; (4) a BEC phase, with long-range coherence of 
spin quasiparticles homogeneously spread over a percolating 
region. In addition to these phases, the Mott insulator (not 
shown), is characterized by absence of excess spin quasiparti- 
cles. The pictures shown here refer to the regime h ~ h c i- An 
analogous succession of phases, mirrored along the magnetic 
field axis, is observed for h ~ h C 2- 

of the critical temperature which obey the same scaling 
law, i c ~ -h)+°. 

In the presence of site dilution the phase diagram 
changes substantially. The results of an extensive quan- 
tum Monte Carlo study (20] for a dilution of x =15% are 
shown in Fig. [T] The most important effect of disorder 
on the phase diagram is that of shifting the critical fields 
for condensation and shrinking the condensation region 
with respect to the clean system, h c \ w 1.94(1) > nfi 
and h C 2 ~ 8.52(2) < ■ This effect separates the onset 
of condensation from the closing of the spin gap, which 
still occurs at and . In fact, the spin gap of 
the doped system is the same as for the clean system, 
given that the minimum spin gap in the doped system is 
associated with local excitations in rare regions that lo- 
cally approximate the clean system, and whose size can 
be arbitrarily large. Hence disorder "splits" an ordinary 
quantum critical point into two, i.e. a spectral transition 
point for the closing of the gap, and a true critical point 
for the onset of order. Inbetween these two critical points 
an extended Bose-glass phase opens up, characterized by 
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the joint absence of a spin gap and of superfluidity. In 
the magnetic language, this represents a gapless quantum 
paramagnetic phase. In the clean case of the Hamilto- 
nian given in Eq. ([l]) the T — quantum-critical point 
is Gaussian, with mean-field critical exponents and dy- 
namical critical exponent z — 2 [5 r In particular, the 
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spin gap closes as S ~ \h— ^1(2) i — i" — '"ci(2) I 
v = 1/2 is the mean-field correlation length exponent. 
In the disordered system the gap still closes in a mcan- 
field-like fashion, but disorder completely discards the 
mean-field picture for the true quantum-critical points 
at h c \ and h C 2i which obey the exponents of a new uni- 
versality class pertinent to the superfluid-to-Bose-glass 
transition. The prediction of Ref. [5] for the dynamical 
critical exponent, z = d — 3, is well verified by our nu- 
merical data. Moreover, we determine the correlation 
length exponents as v = 0.7(1), consistent with Ref. [2"51 
and the order-parameter exponent as f3 = 0.9(1) |22J. 

The Bose-glass state is characterized by Anderson lo- 
calization [19 of the excess spin triplets (or triplet holes) 
injected into the Mott insulating state: this leads to a 
magnetization profile characterized by magnetized "pud- 
dles" , located in the rare regions devoid of lattice vacan- 
cies, as shown by direct inspection of the numerical data 
(see Fig. |2|. To characterize the density profile of the 
localized spin triplets (triplet holes) we investigate the 
behavior of the magnetic participation ratio (mPR) 



mPR = 
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representing the effective fraction of the system which 
is magnetized, Here (...) disorder represents the disorder 
average. Fig. [3] show the mPR as a function of temper- 
ature for various magnetic fields. A percolating magne- 
tization profile is characterized by a participation ratio 
which exceeds a fraction p*(x) = p sc /(l — x) = 0.366, 
corresponding to the percolation threshold for the sim- 
ple cubic lattice p sc = 0.312 [21] renormalized by the 
dilution. Below this threshold, the magnetization profile 
is composed of spatially disconnected puddles, which are 
phase-incoherent, resulting in the absence of long-range 
order of the transverse spin components. This remains 
true for a significant portion of the finite-temperature re- 
gion above the zero-temperature Bose glass, despite the 
strong thermal activation of spin triplets/holes. We term 
this temperature region as thermal Bose glass. In such a 
region the magnetization m(t, h) = J^i (£f)/[(l ~ x )^\ 
can be well reproduced by a model of disconnected clus- 
ters with sizes obeying the statistics of rare clean regions 
in a site-diluted cubic lattice [2"2"| . 

Increasing the temperature above the characteristic 
value t m PR=p* at which the magnetic participation ra- 
tio reaches the condition mPR = p* , one enters into a 
thermal percolation crossover region, marking the upper 
bound of the thermal Bose glass. Above this region the 
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FIG. 3: Magnetic participation ratio (see definition in the 
text) as a function of temperature for different magnetic fields. 
The two thresholds mPR= 0.366 and mPR= 0.5 are indicated 
by dashed lines. 

spin triplets / holes, albeit incoherent, are delocalized 
over the entire lattice. Thermal percolation is marked by 
a substantial change in the statistics of the local mag- 
netization (Sf), corresponding to the local excess pop- 
ulation of triplets / holes. Indeed one can easily show 
that mPR = [1 + (er/m) 2 ] -1 , where m is the average of 
the (Sf) distribution (namely the total magnetization) 
and a its standard deviation. In particular the link be- 
tween mPR and the local magnetization distribution sug- 
gests a second criterion for thermal percolation. Consid- 
ering the ratio R = m/a, we take R = 1 (corresponding 
to mPR= 1/2) as the crossover value above which the 
fluctuations are dominated by the mean. This crossover 
marks the passage from a highly inhomogeneous magneti- 
zation profile to a homogeneous one. Hence in the follow- 
ing we adopt the two criteria (mPR= p* and mPR= 1/2) 
to mark the thermal percolation crossover. Most impor- 
tantly, in absence of spontaneous ordering transverse to 
the field, the (Sf) distribution is probed by the lineshape 
of nuclear magnetic resonance |25j . which means that 
thermal percolation can be experimentally detected. 

Given the above analysis, it is to be expected that the 
transition to triplet/hole condensation, induced by the 
applied field, will be very different when coming from 
the thermal Bose glass or when coming form the normal 
Bose gas. In the former case, indeed, the triplet/hole gas 
undergoes percolation of quantum coherence, and hence 
condensation, via tunneling of the triplet/hole bosons be- 
tween the magnetization puddles; this clearly represents 
a finite-temperature quantum percolation phenomenon. 
In the latter case, in contrast, condensation occurs on 
the magnetized backbone of the diluted lattice, akin to 
field-induced condensation in a clean system. 

As shown in Fig. [4] this crossover from low- 
temperature quantum percolation to more conventional 
condensation on a random network has a striking signa- 
ture in the central feature of magnetic BEC, namely the 
scaling of the critical temperature close to the lower crit- 
ical field. The Bose-glass nature of the disordered phase 



4 




order, and more generally for theories of the fundamental 
interplay between quantum critical fluctuations and geo- 
metric randomness. 
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FIG. 4: Percolation phase diagram of site-diluted coupled 
5 = 1 chains in a magnetic field. In the inset: log-log plot of 
t c vs \h — h c i | . A completely analogous behavior is observed 
close to the upper critical field h c i- 



at low temperature leads to the breakdown of the scaling 
prediction stemming from mean-field theory of BEC in 
a dilute Bose gas. Numerically we find that the onset 
of T c follows t c ex \h — h c i(2)\ with <f> — 1.16(5), in vio- 
lation of the mean-field prediction [22] ■ But, even more 
remarkably, above a characteristic temperature range the 
t c (h) curve shows a clear crossover to a different curva- 
ture, consistent this time with the mean-field exponent 
4>o = 2/3. We observe that this crossover corresponds to 
the dividing region between the two different condensa- 
tion regimes, namely from a thermal Bose glass and from 
a normal Bose gas. Indeed, above the crossover region 
spatial percolation of the magnetization profile occurs 
already before condensation, and hence condensation ac- 
quires a conventional character, while disorder only plays 
the role of effectively renormalizing the Hamiltonian pa- 
rameters of the clean system. It is important to stress 
that the disorder-induced crossover in the scaling behav- 
ior is not accompanied by a change in the universality 
class of the transition: indeed, according to the Harris 
criterion |26j disorder is not pertinent for the 3D XY 
condensation transition of the clean system, for which 
v PS 0.67 > 2/D = 2/3 [27]. 

In conclusion, we have shown that a realistic disor- 
dered quantum spin model, corresponding to the mag- 
netic Hamiltonian of DTN [IT] under non-magnetic dop- 
ing, can elucidate the universal features of the physics 
of strongly correlated bosons on a disordered lattice in 
the grand-canonical ensemble. A variety of measure- 
ment tools, such as magnetometry and NMR, can reveal 
the central features of Anderson localization of triplet 
bosons, leading to inhomogeneous magnetization pro- 
files, and to thermal and/or quantum percolation thereof. 
These results clearly demonstrate the potential of model 
quantum magnets as a valuable testbed for theories of 
complex Bose systems with strong interactions and dis- 
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